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We  consider  quasilinear  hyperbolic  partial  differential 
equations  modeling  ideal  gas  flow  under  various  physical  effects. 
When  these  effects  are  represented  as  Lipschitz  continuous 
functions  of  the  states,  solutions  to  the  initial  value  problem 
are  shown  to  exist  globally  in  time.  Our  analysis  is  based  on 
the  random  choice  method  which  generalizes  the  Glimm  scheme  for 
hyperbolic  conservation  laws. 

When  the  effects  are  strongly  dissipative  then  the  flow 
decays  exponentially  to  a  constant  state  as  time  tends  to 
infinity. 
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SIGNIFICANCE  AND  EXPLANATION 


We  study  the  initial  value  problem  for  quasilinear  hyper¬ 
bolic  system  of  the  following  form: 

(-) 

=  U  (u,v,x,t) 
lv  -  -f^  =  V  (u,v,x,t)  . 

ot  a  X 

The  system  models  ideal  gas  flow  where  U  and  V  represent 
physical  effects  such  as  damping  and  external  moving  force. 
Suppose  that  the  initial  data  (u(x, 0) , v(x, 0) )  have  bounded 
total  variation,  and  U  and  V  are  Lipschitz  continuous 
functions  of  u  and  v.  Then  the  initial  value  problem  has 
an  admissible  weak  solution.  When  U  and  V  represent  strong 
dissipative  effects,  then  the  solution  decays  to  a  constant 
state  as  t  tends  to  infinity.  Our  analysis  is  based  on  the 
random  choice  method  which  generalizes  the  Glimm  scheme  for 
conservation  laws. 

Our  assumptions  are  general  enough  to  include  several 
important  physical  effects.  Moreover,  the  random  choice 
method  can  be  applied  to  more  general  quasilinear  hyperbolic 
systems  as  an  effective  numerical  scheme. 


The  responsibility  for  the  wording  and  views  expressed  in  this 
descriptive  summary  lies  with  MIC,  and  not  with  the  author  of 
this  report. 
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AN  INHOMOGENEOUS  QUASILINEAR  HYPERBOLIC  SYSTEM 
Ching-hua  Wang 


INTRODUCTION 


We  consider  the  initial  value  problem  for  the  inhomogeneous 


hyperbolic  system 


-  O  (u,v,x,t) 
.  “  v 


(0.1) 


(u,v)  (x,o)  *  (UQ(x),  Vq(x)) 


(0.2) 


Here  u  (x)  and  v  (x)  are  bounded  functions  with  bounded 
o  o 

total  variation,  vq (x)  has  a  positive  lower  bound.  It  is 
assumed  that  there  exists  a  constant  M  >  0  such  that  U  and  V 
do  not  depend  on  x  for  |x|  >  M  .  U  and  V  are  smooth 
functions  of  u,v#  x  and  t  . 

By  the  Riemann  invariants 


r  *  u  +  log  v  s  *»  u  •*  log  v 


(0.3) 


the  system  (0.1)  can  be  diaglized  as 


Sr  _  1  Sr 
7t  ”  v  Sx 

Ss  ,1  Ss 
Tt  +  v 


R(r,s,x,t) 


S(r,s,x,t) 


(0.4) 


where  R(r,s,x,t)«U  +  -  ,  S(r,s,x,t)  -  V  -  ~  ;  (U,V)  =  (U,V)  (£~,exp (£j£.),x, t) 

Throughout  this  paper  we  also  assume  that  there  exist  two 
constants  K^»  K2  such  that 
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iRfrj^s^Xft)  -  R(r2,s2,x,t)  |  Kx  ( |  rx-r2 1  +  Isj^-Sjl) 
|S(r1,s1,x,t)  -  S  (r2,s2,x,t)  |  <  K1  { | r^— r2 1  +  1  si“s2 1 ) 


(0.5) 


|£|  < 

dX1  - 


fl 


(0.6) 


When  (U,V)  =  (0,0),  system  (0.1)  reduces  to  the  following 


system 


9u 

3t 


0 


3v 

Tt 


(0.7) 


The  initial  boundary  value  problem  of  (0.7)  has  been  solved  by 
Nishida  T.  [6] .  When  U  and  V  do  not  depend  on  x  and  t,  the 
global  solution  of  (0.1  and  (0.2)  has  been  constructed  by  Ying 
Lung- an  and  Wang  Ching-hua  [7] .  For  general  inhomogeneous  system 


9u 

3t 


+ 


3f  (u)_ 
3x 


g(x,t,u) 


(0.8) 


Ying  Lung-an  and  Wang  Ching-hua  18)  have  proved  the  existence  of 

the  global  solution  of  initial  value  problem  of  system  (0.8)  for 

_Kt 

g(x,t,u)  *  g(u)  e  .  When  g  does  not  depend  on  t,  Liu  Tai-Ping  [5] 
has  constructed  global  solutions  for  system  (0.8)  and  studied  their 
asymptotic  behavior. 

In  general  (0.1)  and  (0.2)  does  not  possess  smooth  solutions, 
and  we  look  for  weak  solutions  in  0^t<T,  i.e.,  solutions  satisfying 


3. 


+0B 

II  (u<J>t  +  7  <t>x  +  u4>  )  dxdt  +  I  Uq(x)  <Mx#0)  dx 


0<t<T 


II  (vV»t  -  u\|>x  +  Vifi)dxdt  +  |  Vg(xH(x,0)  dx  =  0 


0<t<T 


(0.9) 


(0.10) 


for  any  smooth  functions  $  and  with  compact  support  in  0<_t<T. 

The  purpose  of  this  paper  is  to  prove  the  existence  of  the 
weak  solution  for  (0.1)  and  (0.2)  and  study  their  asymptotic 
behavior  as  the  time  variable  t  tends  to  infinity. 
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1.  Difference  Scheme 


We  now  describe  the  difference  scheme  we  use  for  solving 
(0.1)  and  (0.2),  which  is  the  generalization  of  one  in  [2]  and 
[7J.  Randomly  choose  an  equidistributed  sequence  a={ai>  in 
(-1,1)  and  mesh  length  Ax«=.l  At«h  satisfying  the  Courant- 
Friedrich-Lewy  condition 

l  1 

^  -  const  >  - 


for  all  v  under  consideration. 


Let  (u,v) (x,t,  m+1,  k,  u£,  v£,  ur,vr))  denote  the  solution 
of  the  system  (0.7)  at  the  point  ( (m+1)  l,  kh)  with  the 
following  initial  data 


(u,v) (x,kh)  =. 

K 


v^)  x  <(ro+l)  l 
vr)  x  >(m+l)  l 


(1.1) 


We  take 

(u,v)  (x, t ,  m+1,  k,  u,  ,  v4  ,  ur,vr)=  (u+U(u,v, (m+1) i,  kh) (t-kh) , 


v  exp  (t-kh))  (1.2) 

v 

as  the  approximation  solution  of  the  system  (0.1)  with  the  initial 
data  (1.1).  In  fact,  using  the  Riemann  invariants,  we  can  rewrite 
(1.2)  in  the  following  form 

(r,s)  ■  (r+R(r,s,  (m+1)  £.,kh)  (t-kh)  ,2+S  (r,S,  (m+1)  i,kh)  (t-kh) )  (1.3) 


5. 


The  construction  of  the  difference  solution  (u,v) (x,t) 
proceeds  as  follows: 

For  each  mesh  length  l  ,  let  (u.  (x)  ,  v.  (x) )  be 

x  #o  X  /  o 

defined  by  the  equation 

■  tuoc-»)  ,(vo(-»n ,  x  <  -(2i|-]2  +  m 

(u4fo  v£,o)(x)=  '  C^<W).v0(W>>,  for | x|  <(2Ij]2+1)£  and 

Cm-1) l <x<(m+l) l,  m=even 

(u0(+»),v0  (+“)),  x  >  (  (2  I 2  +  1)  £  (1.4) 

Then  we  set 

(u,v)  (x,t)  *  (u,v)  (x,t,  m+l,o,  (u£  o,vA^(mJl)  ,  (u^,^  )  (m+2)  S.) ) 

for  mi  <  x  <  (m+2)£,  0  <  t  <  h,  m  =  even  (1.5) 

Inductive  ly  suppose  that  (u,v) (x,t)  has  been  defined 
for  0  <  t  <  kh,  we  set 

(u,v)  (x,t)  =  (u,v) (x,t,  ro+l,k,  (u,v)  ( (m+akU  ,  kh-0)  (u,v)  ( (m+2+akU  ,kh-0)  ) 

for  mJt  <x  <  (m+2)i,  ,  kh  <  t  <  (k+l)h,  m+k  =  even  (1.6)_ 

It  is  obvious  that  (u,v) (x,t)  depends  on  the  mesh  length  h  and 
sequence  a  but  we  omit  the  subscript. 


2.  Bounds  for  the  Difference  Equation 

Because  we  use  (u,v)  defined  in  (1.2)  as  the  building 

block  to  build  the  difference  solution  of  (0.1)  and  (0.2),  we 

have  to  do  some  estimates  on  it. 

First  we  consider  (u,fr)  (x,t,  (m+1) ,  k,  UA'V£»  ur'vr* 

defined  in  Section  1.  It  is  well  known  that  there  are  three 

constant  region  in  the  half  plane  t  >  kh,  in  which  (u,v)  equals 

to  (uA#vJl),  (Vvm)  and  (ur,vr)  respectively.  (um'vm)  is 

connected  to  (u0,v.)  by  a  1-wave  on  the  left  and  (u  ,v  )  is 

x  x  r  r 

connected  to  (u.v  )  by  a  2-wave  on  the  left.  We  define 
mm 

(r.,s.),  (r  ,s  )  and  (r  ,s  )  by  means  of  (0.3)  and  define 
xx  mm  r  r 

(2.1) 
(2.2) 

(2.3) 


(2.4) 


Ar  =  r  -  r  As  =  s„  -  sm 

m  x  r  m 


P((ut,vA),  (ur,vr))  ■  -  min (o, Ar)  -  min  (o,As) 


It  has  been  proved  [6]  that 


P(  (u^Vj) ,  (Uj,Vj) )  <  P((ui,vi)r  (uk»vk))+P((uk,vk),  (UjfV^)) 


where  u^,u^,Uj  are  any  constants,  vi,vK,vj  are  any  positive 
constants 


Let 


<r  *  rr  "  ri  '  Ss  *  sr  ‘  S1 


we  have 


6r.  «=  u  +  U  (t-kh)  +  log  (v  exp  -^(t-kh) ) - (u.+U# (t-kh) +log (v.exp_^(t-kh) 
trr  rv  *  *  * 

r  vi 

■  <rr"rj,)  +  CRr-  R£)  (t-kh) 

vr  V£ 

«s.  =  u  +  U (t-kh) -log (v,  exp  (t-kh) )- (up+Up (t-kh) -log (vfexp— -(t-kh) ) 

trr  r  ^  v  £ 

-  (sr  -  Sjl)  +  (Sr  -  8t)  (t-kh) 

where  U  ,V  ,R  ,S,.  evaluated  at  (u,v  f  (m+1)  trkh)  and  Up,  Vp  , 
rrrr  *  r  a#*# 

Rl'  S£  evaluated  at  (u^v^,  (^-+1)  kh)  . 

Therefore  by  (0.5)  we  have 
Lemma  2.1 

1 6r.  -  6r|  <  K.(|6r|  +  |6s|)(t-nh) 

t  "  1  (2.5) 

|6s  -  6s |  <  Kx ( 1 6r |  +  J  6s | ) (t-nh) 

In  the  same  way  as  Lemma  3  in  [7] ,  we  have 

Lemma  2.2  If  (ur»vr^  can  be  connected  to  (u^v^) 
only  one  wave,  u*,  u*  are  any  constants  and  v*,  v*  are  any 
positive  constants  then 

P((u*,v*),  (U*,v*  ))  <  P((uifvt),(ur,vr))  +  3(|6r*-6r|+|6s*-6s| )  (2.6) 

where  6r*=r*-r £= (u*+  log  v*)-(uj+  log  v*),6*=(u*-  log  v*)-(uj-  log  v*) . 

In  order  to  estimate  the  variations  of  r(x,t)  and  s(x,t), 
we  define  some  functionals. 

For  integers  k  1  ,  we  define 

F(kh+o)  «  l  P((u,v)  ( (m+a.  )  l ,  kh-o) ,  (u,v)  ( (m+2+a.)  i,kh-o) )  , 
m+k^even  * 


and  for  integer  k  >_  0 

F((k+l)h-0)  =  l  (P( (u,v) ( (ro-1  +  a,  ,) £, (k+l)h-0) , 
m+k-even 

(u,v)  (m£.-0,  (k+1)  h-0) )  + 

P ( (u,v) (m£-0 , (k+1) h-0) , (u,v) (m£  +0  ,  ( (k+1)  h-0) )  + 

P  ( (u,  v)  (m£+0 ,  (k+1)  h-0),(u,v)  ( (m+l+ak+1)  l  ,  (k+1)  h-0) )  } 

4-°°  +a> 

Lemma  2.3.  For  any  given  T  >  0,  V  u(*,t)  and  V  v(*,t)  the  total 
—————  —  00  —00 

variation  of  the  difference  solution  (u(x,t),  v(x,t)),  are 

bounded  uniformly  for  the  mesh  length  h,  sequence  a  and  t  <  T, 

v(x,t)  has  positive  lower  bound  uniformly  for  the  same  parameters 

provided  k<  ^  . 

2K1 

Proof  .  Since  U(u,v,x,t),  V(u,v,x,t)  do  not  depend  on  x  for 
I x|  >  M  and  (0.5)  there  exists  a  constant  M  depending  only  on 
T  and  uQ(±°°),  vQ(±oo)  such  that  |s(±»,t)|  <  M,  |r(±«,t)|  <M  for 
0  <  t  <  T#  provided  h  is  small  enough* 

It  follows  from  that  U(u,v,x,t) ,  V(u,v,x,t)  do  not  depend  on 
x  for  |  x|  > M  again  and  (0.6)  that 

l  P( (u,v)  (m£-0, (k+1 ) h-0 ) , (u,v) (m£+0,  (k+l)k-0) )  < 

m+k=even  ~ 

h  l  (lR(  (r,s)  (m+a.  )£  ,kh-0)  ,  (m-l)i  ,kh)  -  R(  (r,s)  (m+a.)£  ,kh-0) 
m+k=even  K 

( ra+l)£  ,kh)  |  +  |s ( (r,s) ( (m+ak)£  ,k  -0) , (m-1) ,kh) -  S ( (r,s)  (m+ak)£  ,kh-0) 

(m+l)2.  ,kh)|}  _<  4K2  Mh  . 
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+  00  +00 

It  is  obvious  that  V  r(*,nh+0)  +  V  s(*,nh+0)^  4F(nh+0)  +  4M. 

—  00  — oo 

Then  from  Lemma  2.1,  Lemma  2.2,  we  obtain 


F  l  (n+1) h-0)  <  F  (nh+0)  +  eKjhCV  r(*,nh+0)  +  V  s(*,nh+0))  +  4K2Mh 


_<  F(nh+0)  (l+24K1h)  +  (24K1Mh+4K2Mh) 


Because  F((n+l)h+Q)  <  F( (n+1) h-0)  we  have 


F  ( (n+1)  h+0)  <  F  (nh+0 )  (l+24K]Lh)  +  24K1Mh+4K2Mh) 

Let  F (0)  =  l  P( (uQ,v0) (mH) , (uQ, ( (m+2) l) ) ,  then  we  obtain 
m=even  u 

+  00  +GO 

F(0)  <  V  rQ(*  )  +  VsQ(-),  where  rQ(x)  =  uQ  (x)  +  log  vQ  (x) , 

—00  —00 


sft(x)  =  u-(x)  -  log  vn(x). 


After  simple  computation  we  get 

K,M  24K.T 

F(  (n+1)  h+0)  <  (F  (0)  +  ft  +  )  e  A 

+  »  +  ®  _  K.M  24K.  T 

<  (V  rQ  ( • )  +  V  3.(0  +M  +  )  e  1  for  (n+1 )  h  <  T 

“  OO  —  OO  2. 

and  thus 

®  +  eo  +00  +O0  ^  K2M  24K.T 

V  r  (•  ,nh+0)  +  V  s(*,nh+0)  <  4  (V  rQ(*)  +  V  sQ(*)  +  M  +  )e  1  + 

■00  —  00  —  oo  —oo  1 


+  4T1  for  nh  <  T 
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It  follows  from  that  Uq(x)  and  Tq(x)  are  bounded  functions 

with  bounded  variation  and  vq (x)  has  a  positive  bound.  Lemma  2.1 

and  h  <  that 

2K1 


V  r(.,t)  +  V  s(.,t)  <  Q 


for  t  <  T 


(2.7) 


+oo  +oo  j {  jgj  24K  T  K 

where  Q  =  8  (  V  rft(«)  +  V  sft(«)  +M+  )e  1  +  ^  M  +  8U  . 

-co  0  _»  0  6K1  K1 


Therefore 


|r(x,t)  |  <  2M  +  Q,  |  s  (x,t)  |  <  2M  +  Q,  |  u(x,  t)  |  <  2M  +  j  for  t<T  (2.8) 


-  (2M+  2)  2M+  % 

e  <  v(x,t)  <  e 


+00 

v  u(  *,t)  <  Q/2 


+«  ,  2M+  2- 

W(.,t)  <  jO  e  z 
—00 


for  t  <  T  (2.9) 


for  t  <  T  (2.10) 


for  t  <  T  (2.11) 


Q.E.D. 


Lemma  2.4.  For  any  given  X  > 0,  there  exists  constant  D  depending 
only  on  U,V,  uQ(x),  Vg(x),  T  and  X  such  that 
X 

I(X)  «  j  (|u(x,t2)  -  u(x,t1)  I  +  |v(x,t2)  -  vfx,^)  |)  dx  < 


<  D((t2~t1)  +  h)  ,  0<t1,t,<T  (2.12) 


provided  h  is  small  enough. 
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Proof:  For  definiteness  suppose  and  there  is  a  pair 
of  integers  n,k  such  that  nh  t1  <  (n+l)h  <...<  (n+k)h<  t2 
<  (n+k+l)h.  It  is  obvious  that 


I(X)  <  I X(X)  +  I2(X) 


where 


I^X) 


I2(X) 


k+1  X 

l  f  ( |u(x, (n+i)h+0)  -  u(x, (n+i)h-0|  + 

i«°  Jx 

+  |  v(x,  (n+i)h+0)  -v(x,  (n+i)h-0)|  )  dx, 

k  X 

l  (  (|  u(x,  (n+i+l)h-0)  -  u(x,  (n+i) h+0)  |  + 
i=0  > 

- X * 

+  |  v(x,  (n+i+1)  h-0)-v(x ,  (n+i)  h+0)|)  dx 


X 

+  J  (lutx/tj^)  -  u (x,nh+0)  I  +lv  (x,t^)  -  v (x,nh+0)  |  )dx  + 

-X 

X 

+  |  (|u(x,t2)  -  u(x,  (n+k)h+0)  |  +  |v(x,t2)  -  v (x, (n+k) h+0) | ) dx. 
-X 


It  is  not  difficult  to  show  that 


I^X)  <  21 


*  4*00  +00 

sup  _v  u(*,t)  +  sup_Vv(*,t) 

{  0<t<T  0<t<T 


For  x  e ( (m-l)i,, (m+1) l) ,  m+n+i+1  =  even.  Using  (1.2)  and  (2.9) 
we  have 
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(TO+Dl  £ 

|u(x,  (n+i+l)h-0)  -  u(x,  (n+i)h+0)  |  <,  V  u  (x,  (n+i)  h+0)  +  |u  |h. 

(m-l)4 

(xn+1)  i  * 

[ y(x, (n+i+1) h-0)  -  v(x, (n+i) h+0)  |  <  V  r (x, (n+i) h+0)  +  2  |v  |h 

”(ro-m 

provided  h  is  small  enough,  where 


(U*,V*)  -  ((OfVj(x,  (n+i+1) h-0)  ,  m*,  (n+i)h)  . 

X 

For  I  (|  utx,^)  -  u(x,nh+0)|  +  |v(x,t)  -  v(x,nh+0)|  )dx  and 
*-X 

(|  u(x,t2)  ~  u(x,  (n+k)h+0)|  +  |v(x,t2)  -  v  (x,  (n+k)  h+0)  |  )  dx 
JL 

we  have  similar  estimate. 

Combining  the  above  estimates,  we  obtain 


+  00  +00 

I-  <  (2  sup  V  u(*,t)  +  sup  V  r(*,t)  + 

*  [  0<t<T  -oo  0<t<T 


* 

2  sup  V 
m+n+i+l=even 


Since  (2.8)  (2.9(2.10)  and  (2.11),  we  obtain  (2.12).  Q.E.D. 
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3.  Convergence  of  the  Difference  Solutions 

Because  (u,r) (x,t)  is  not  an  exact  solution  of  (0.1)  in  the 

T 

strip  kh  j<  t  (k+l)h,k  =  0,1,...,  jj-  ,  so  in  order  to  prove  the 
existence  of  the  weak  solution  of  (0.1)  and  (0.2)  we  first  prove 
two  lemmas. 

Lemma  3.1.  If  x'  =  £^t'  is  a  l-schock  wave  of  (u,v) 

(x,t, (m+1) ,k, (uA,v^) , (ur,v  ) )  in  che  domain  {(x,t)s  m  <  x  <(m+2)£.  , 
kh<  t  <  (k+l)h,  m+k  =  even)  then 


I  €1Iu]  -  |<  wt*  (1%-uJ  +  |vm-vj)  (3.1) 

ICjIv]  -  lu]  I  <  wt'(  u^-uj  +  v^-vj)  (3.2) 

provided  h  is  small  enough,  where  x'  ■  x  -  (m+l)&,  t'  =  t  -  kh, 
the  constant  w  depends  only  on  U,V,  Uq(x),  rg(x);  (u^,v&)  *  (u,v) 


( (m+a^)  i,kh-0),(ur,vr)  =  (u,v) (m+2+a^) &,kh-0,  ;  (um/Vm)  is  the 
constant  solution  being  connected  to  by  1-shock  wave  on 

the  left;  J  J  denotes  the  jump  of  the  quantity  in  bracket  across 
the  1 -shock  wave. 

For  x'  **  ^2^'  is  a  2-shock  wave  of  (u,v)  ,  we  have  same  kinct 
of  result. 

Proof:  It  follows  from  that  (u  ,v_)  and  (u_,v_)  satisfy  Rankine- 

—  ’  1  r  r  mm 

Hugoniot  condition  and  (0 . 5)  (2 . 8)  (2 . 11)  that 


14. 


vt* 

in 


5  l«1<VVt'l  *  I  7,  / 


e-r|  dn 


l  Jo 


1  (  v" 

%  Jo 


e  n  dn 


v*t'  v.f  v 

i  U1<VVt'l  +  I  <vt-V  <-  ~2  f  v*  +  c  •“  17  (v7»l  + 

V*  'o 

Vt' 

+  1  v„  /*  e'nd"li  wt’<lum-all  +  lvm'vJ> 

®  V  t’ 


where  (U^/V^)  =  (U,V)  (u^v^,  (m+1)  A,kh)  ,  (Uro,Vm)  =  (U,V)  (m+1)  l ,kh) 

V*  **  V(u4,v*,  (m+1)  A,kh)  ,  V  =  V(ujt,vjn,  (m+1)  Jl,kh)  ; 

v*  ■=  8um  +  (l-e)vA,  O<0<1  .  The  proof  of  (3.2)  is  similar.  Q.E.D. 

Lemma  3.2.  If  (u,v)  is  a  1-rarefaction  wave  in  the  subdomain 
*  {(x,t)  ;  £q  <  -T  <£^f  0<  t'<  h}  of  the  domain  fi  =  {(x,t), 

mi  <x  <(m+2)£,kh  <  t  <  (k+l)h,  m+k=even)  ;  then 


I  ( (u4>t  +  v  ♦x  +  u>dxdt  +  |  ufrds  -  ±  <J>dt| 


3  2  C»t* 

<  R.h  +  R,h  (  V  u (kh+0)  +  V  v(.  ,kh+0)) 
“  A  *  r»  t*  £•  t' 


(3.3) 


1 1  (v  Xt“  uxx+V)dxdt  +  |  vxdx  -  ux<Jt| 


3  2  S.t* 

<  R.h  +  R,hz(  V  u( • ,kh+0)  +  *V  v(*,kh+0) 

1  *  C'«f  ('At' 


(3.4) 


15. 


provided  h  is  small  enough. 

Where  x*  ■  x  -  (m+l)H,  t=t-kh;  (u,v)  =  (u,^)  (x,t,  (m+1)  ,k,  (u, v) 
(ro+ak)  «-,kh-0)  ,  (u,v)  ( (m+2+a^)  l,  kh-0) )  ,  <p  and  x  are  smooth 
functions  with  compact  support  in  0^t<T  ;  R^  and  Rj  are  corns  tan  ts 
depending  only  on  U,V,Uq(x),  (x)  ,  <P  and  x  »(U#V)  =  (U(u,v,x,t)  , 

V  (u,v,x,  t) ) . 

For  (u,v)  is  a  2-rarefaction  wave  .in  the  subdomain 

« 

ft  2  “  { (x,t)  ;  £'  <  ?2 '  ^<t<h  }  of  the  domain  ft  ,  we  have  same  kind 

of  result. 

Proof:  (u,v)  ,  as  a  1-rarefaction  wave  in  ft^,  satisfies  (0.7)  in 

in  classic  sense.  Because  of  (0.5)  and  (2.9),  U  and  V,  as  the 
functions  of  u  and  v,  satisfy  Lifshitz  condition  with  the  Lifshitz 
cons tern t  depending  only  on  ,  M  and  Q. 

It  follows  from  (1.2) (1.6)  that 

3u_.3u  3U3u.3U3v..,  /-j  c\ 

3x  "  3x  +  3~3x  ^~3x,t 


3u 

3t 


3U 

3u 


3u  3U 
3v 


3v 

3t 


)t'  +  U 


3v 

•§3T 


vt* 

v 

e  (v 


(3.6) 


+  v 


3V 

3V 


ii  _  V 
3x 


3v 

3x 


)  (3.7) 


.  iViMiiiiaiiifiiifiii  in  i 


|v=|vev  +  (v,-lf  ev  (vf  +  +ve 


where  (U,V)  =  (U,V) (u,v, (m+1) A,kh)  . 

8U  s4> 

Let  -  U,  then  we  obtain 


(3.8) 


(Vt'  ) 


CX*’) 


e  +  (v)  x  t‘  e 


U  =  (1-e 


C  vt'  1 

lT  1 


*>'x  *• 


-3  _  * 


fv  4 

\  du 


3u  ~  3V  3v  3v  \ 

"3x  +  v  "5x  v  "5x  Jr 


Using  (0.5)  and  (0.6)  for  estimating  U  -  U,  then  we  get 


ff  (u<|>.  +  ^  4>  +  U$)dxdt  +  f  u<j>dx  -  ^  4»dt 

})  t  v  x  j  3n^  v 

fti  ?lt’ 

<  R.h3  +  R_hY  V  u(*kh+o)  +  V  v(*,kh+0)\ 

-  1  2  \«f  qf  ; 


The  proof  of  (3.4)  is  similar. 


=  ||  (x ,  t )  4»  (x ,  t)  dxdt  < 


Q.E.D. 


We  use  the  difference  scheme  introduced  in  section  2  and  take 

2ff+0  m 

the  ratio  of  the  mesh  lingths  r— ~  “  e  w  and  h.  *  — .  .  When 

hi  1  21 

i  is  large  enough,  the  difference  solution  (u,v) (x,t)  can  be  defined 
in  0<t<T. 
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The  difference  solution  (u,v(x,t)  is  discontinuous  at  the 
segment  {(x,t):  x  ■  mi.^,  kh^  <  t  <  (k+l)h^,  m+k*even)  .  But  it 
follows  from  (0.6)/  (2.8)  and  (2.11)  that 

(k+l)hi  (k+1)!^ 

if  (u(m  .+0,t)  -  u(ro£  .-0/t)  )x  dt|  <  f  (t-kh.)|  U  ( (u,v)  ( (m+a.  )  l 

‘  kh±  1  1  ~  Jkh±  x 

(m+1)  i^kh^  -  U(  (u,v)  ( (m+ak)  (ra-1)  fc^khj  x  dt  <  E  h|  (3.9) 


provided  i  is  large  enough,  where  E  is  a  constant  depending 
only  on  K2,  M  and  Q. 

With  Lemma  3.1,  lemma  3.2  and  (3. 9), (3. 10)  in  mind  by  the 
way  similar  to  these  in  the  section  4  in  [7]  we  obtain  the 
following  existence  theorem. 


Theorem  3.1.  Suppose  that  Uq(x)  are  bounded  functions  with 
bounded  variation,  Vq(x)  has  a  positive  lower  bound,  U(u,v,x,t), 
V(u,v,x,t)  are  smooth  functions  satisfying  (0.5)  (0.6)  and  do  not 
depend  on  x  as  |x|>M,  where  M  is  any  given  positive  number, 
then  for  any  given  T>0  the  initial  value  problem  (0.1)  and  (0.2) 
has  a  weak  solution  in  0<  t<T. 
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4 .  Asymptotic  behavior 


Now  we  consider  the  asymptotic  behavior  of  the  solution  of 
(0.1)  and  (0.2)  under  the  following  additional  restrictions  on 
U  and  V: 


R  depends  only  on  r, 

and 

R'  (r)  -  a; 

(4 

.1) 

S  depends  only  on  s. 

and 

S' (s)  <  -  a 

(4 

.2) 

where  a  is  a  positive  number. 

Lemma  4.1.  Suppose  that  R(r)  satisfies  (4.1),  then  the 
initial  value  problem 


(  It  -  R<r> 


r (0)  =  r. 


(4.3) 

(4.4) 


possesses  a  unique  solution  r(t,rg)  satisfying 


|r(t,r0)-  r |  <  |rQ-r|  e_ott 


for  t  >  0 


(4.5) 


where  r  is  the  unique  solution  of  R(r)  =  0;  and  for 
f 


rh(t'r0)= 


r0+R(rQ) t  ,  0  <  t  <  h 

rh(h,rQ)  +  R(rK (h,rft) ) (t-h)  ,  h  <  t  <  2h 


h  '  0 


,  the  Euler's 


x.  (nh,rn)  +  R(nh,rft)  (t-nh) ,  nh<  t  <  (n+l)k 


difference  solutions  of  (4.3)  (4.4),  the  following  inequality  holds 
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rhtt,r0*  “  rl  1  lr"r0*  e 


-at 


for  t  >  0 


(4.6) 


provided  h  <  —  . 


Proof.  Since  the  assumptions  (0.5),  (4. 3)  and  (4.4)  possesses  a 
unique  solution  r(t,rQ)  for  all  t  >  0.  By  (4.1)  we  have 


d(r (t,rQ) -  r) 
dt 


=  R(r (t,rn) )  =  -a (r (t,rrt) -  r) 


and  therefore  (4.5). 

Now  we  turn  to  prove  (4.6). 

Without  loss  of  generality  we  assume  rg  ^  r,  because  of 
rh(t,r)  =  r(t,r) . 

For  definiteness  suppose  rQ  >  r.  It  follows  from  (0.5)  (4.1) 
that  for  any  r  >  r  and  tQ  >  0  the  initial  value  problem 
dr 


3t 

rl 


R(r) 
=  r 


has  a  unique  solution  r(t,r,tQ)  such  that 


t*»t. 


d  r (t,r, tg) 

dt2 


>  0;  therefore 


r  <  r^tt^g)  <  r(t,rQ) 


for  t  >  0 


(4.7) 


provided  h  <  —  * 


Siroilary  for  rn<  r  we  have 


r(t,rQ) <  rh(t,rQ) <  r 
By  (4.5) (4.7) (4.8)  we  obtain  (4.6). 


for  t  >  0 
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Lemma  4.2.  Under  the  assumption  (4.1)  and  (4.2)  we  have 

F(k+1) h-0)  <  (l-(l-c)  cth)F(kh+0)  ,  for  k  >  0  (4.9) 

provided  h  <  min(-  ^  ) ,  where  c  is  a  constant  depending  only 

+«®  +»  i 

on  V  u^(0,  _V  Vq(*)  ,  the  positive  lower  bound  of  vQ(x)  and 

the  equation  (0.7). 

Proof:  Now  we  consider  (u,v) (x, t, (m+1) ,k,u^,v^,ur,vr)  defined 

in  section  1  again  and  first  remind  some  results  in  [2]  and  [7] . 

2 

There  is  a  c  function  f(x)  satisfying 


f  (x)  =  constant  for  x  _>  0 

0<  f'(x)<  1,  f"(x)  <0,  for  x  <  0  (4.10) 

lim  f ' (x)  =  1 

X-*- 


such  that  is  connected  to  (u  ,v  )  by  1-wave  if  and  only 

if 

sm  “  s*  -  f^m  “  (4*U) 

where  r  -  r„  >  0  means  that  the  1-wave  is  a  1-rarefaction  wave, 
in  x,  — 

and  r_  -  r „  <0  means  that  the  1-wave  is  a  1-shock  wave;  and 
m  x, 

(u  ,v  )  is  connected  to  (u.,v0)  be  2-wave  on  the  left  if  and  only 

XT  IT  X  X 

if 


v  -  r  =  f(s  -sj 
r  m  r  m 


(4.12) 


>  0  (s„  -  <  0)  means  that  the  2-wave  is  a  2-rarefaction 

r  in  —  r  m 

wave  (a  2-shock  wave  respectively) . 
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By  the  notation  of  (2.1)  and  (2.4)  we  have 


6r  ■  Ar  +  f(As)  ,  <5  s  =  As  +  f(Ar) 


(4.13) 


Because  of  (4.10)  we  get 


Ar  =  p(Sr,5s)  ,  As  =  q(6r,5s) 


(4.14) 


where  p,q  e  c  .  It  is  easy  to  know  from  simple  computation 
that 


3p 

~WZ) 


_iP 


1-f ’ (P)f ' (q) 


-f' (q) 

1-f'  (p)f '  (q) 


9(6r) 


IVs) 


-f ' (p) 

l=f '  (p)f  r('q) 


(4.15) 


l=f'  (p)f  *"(q) 


Suppose  that  (u  ,v  )  is  connected  to  (u0,v0)  by  a  1-shock 
in  m  x  x 

wave  on  the  left.  Since  (0.5)  (1.3)  and  (4.1),  then  we  have 


sm -*l  <(sn  +  hs(V>  -<st+hs<8l)> 


ri.-rt  <  <V  hR(rm))  -(rt+h«(rt)) 


<  0  , 


therefore  P(  («t#h.  v^)  ,  C^,h'vm,h,)  <  *<VV  '  ( >  (4.16) 


V(u. ,v. )h 


where  u^h  =  ui+U(ui,vi)h  v^  **  vA  exp(  - — 


)  , 


i  can  be  i,  ro,  or  r. 


1 


22. 


Suppose  that  (u^v^)  is  connected  to  (u^v^)  by  a  1 -rarefaction 
wave  on  the  left.  Since  (0.5) (1.3)  and  (4.1),  then  we  have 


*m  +  hS(sm>  '  al+  S(sl> 


0  <  rB  +  hR(r„)  -  (r,  +  hR(rt) 


therefore 

P<(ut,h'  vi,h>’  (um,h'  vm,h^  ^  '  P“W<  <V  vm>>‘  0  <4‘17) 

The  same  kinds  of  results  as  (4.16)  and  (4.17)  can  be  obtained 
for  a  2-shock  wave  and  a  2-rarefaction  wave  in  (u,v)  respectively. 

By  combining  the  estimates  (4.16) , (4.17)and  (2.3)  we  have 


F((k+l)h-0)  <  F(kh+0) 


for  k  >  0 


(4.18) 


+00  +oo 

therefore  for  JDV  r(»t)  and  DV  s(*,t),  the  decreasing  variation 
—  00  —00 

of  r(x,t)  and  s(x,t)  on  the  interval  (-",  +<*>),  the  following 
inequalities  hold 


+<D 


+00 


+00 


BVrCt)  <F(0)  <  (  Vr.M  +  V  sQ(0)  -  Q-, 

—00  —  OO  —00 

+00  +oo  +00 

DV  s(*,t)<  F(0)<  (  V  rn  ( • ;  +  V  sn(*))  =  Q. 


(4.19) 


It  is  not  difficult  to  know  from  (4.20)  that  the  strength  of 


the  shock  wave  in 

(a,tf) 

is  less  than 

»  i • e . 

|Ar| 

iQi 

for 

1 -shock  wave 

(4.20) 

|As| 

lQi 

for 

2-shock  wave 

(4.21) 

Since  (4.10)  then 

0  <  c  =  f'  (-Q^  <  1  (4.22) 

Now  we  are  ready  to  obtain  the  refinement  of  (4.16)  suppose  that 
(Vv  )  is  connected  to  (u^v^)  by  a  1-shock  wave  on  the  left.  By 
(4.15)  we  have 


vH,h>'(Vh'  Vl.”  '  P<(ul-  V'  (um'V>  * 


/  a?  +/  l-fMp)f*,(q)dn  -  (1'C)  ‘ 

6r  6  s 


(6rh~6r)  + 


+  (6sh-6s).  >  (l-c)h  t(R(rm)  -  R(ra)  +  (S(sm)  -  S  (s^) )  J 


m 


>  -  (l-c)cx  h  Ar 


(4.23) 


where  6r.  **  (r  +  hR(r  ))  -  (r8  +  hR(r0)) 


m 


m 


«sh  -  (s^  +  hS(sm))  -  (s4  +  hS  (s^) ) 

Suppose  that  (ur»  vr)  is  connected  to  (u^/V^)  by  a  2-shock  wave 
on  the  left.  Similarly  we  have 


2 


p“Vh*  Vh>’  (ur,h>  vr,h"  '  Pl(un  'V'^r'V'  i 


>  —  (l~c) a  h  As  (4.24) 

By  combining  the  estimates  (4.23) (4.24) (4.17)  and  (2.3) ,  we  have 


F(k+1)  h-0)  <  (1-  (1-c)  cth)  F(kh+0)  . 


Q.E.D. 


We  turn  now  to  proving  the  existence  of  the  weak  solution 
of  (0.1)  and  (0.2)  and  studying  its  asymptotic  behavior.  We  use 
(u,v) (x,t)  to  denote  the  difference  solutions  of  (0.1)  and  (0.2) 
before  for  convenience,  although  the  difference  solutions  of  (0.1) 
and  (0.2)  depend  on  the  mesh  length  h  and  the  sequence  a.  In 
order  to  distinguish  between  the  difference  solutions  of  (0.1) 
and  (0.2)  and  the  weak  solution  of  (0.1)  and  (0.2),  from  now  on 
we  use  (u*1, v*1)  (x, t)  to  denote  the  difference  solutions  of  (0.1) 
and  (0.2)  and  (u,v) (x,t)  to  denote  the  weak  solution  of  (0.1)  and 
(0.2). 

Theorem  4.1.  Suppose  that  Uq(x),  Vq (x)  are  bounded  functions  with 
bounded  variation  and  do  not  depend  x  as  |x|>  M,  Vq(x)  has  a 
positive  lower  bound.  U,V  are  smooth  function  satisfying  (0.5), (0.6) 
(4.1)  and  (4.2)  then  the  weak  solution  (u,v) (x,t)  of  (0.1)  and  (0.2) 
exists  in  t  >  0  and  satisfies  the  following  decay  laws. 
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V  u(*,t)  <  2QX  e”  at  +  q2  e"at 


(4.25) 


+  00 


V  log  v(*,t)  <  2QX  e"(1_c)at  +  Q2  e"at 


(4.26) 


u(x,t)  - 


r+s 


<  2Ql  e-*1-'101  +  |  Q2  e-at 


(4.27) 


log  v(x,t)  -  ^  <  2QX  e_(1"c)at  +  \  Q2  e~at 


(4.28) 


+  09  +09  __  _ 

where  Q.  =  V  rQ(»)  +  V  sQ(0#  Q2  ■  j(  |  rQ  (-«) -r  |  +  |rQ(+®)-r| 
—00  W  —00 

+  |sq(-«>)-s|  +  |  Sq  (+“>)  -s  | )  ,  r  and  s’  are  the  unique  solution  of 
R(r)  =  0  and  S(s)  =0  respectively,  c  =  f't-Q-.). 


Proof.  If  we  choose  the  ratio  of  the  mesh  lengths 


r*+s* 

r-  *»  const  >  e  2  +  2Q. 


where  r*  =  max (r^ (-“) ,x) ,s*  =  max(SQ (-•) ,s) ,  and  h  ■  2  1  ,  where 
i  is  large  enough  such  that  h  <  min  (^,  ) ,  then  the  difference 

V  W  1 

solution  (u  ,vn) (x,t)  defined  in  section  1  can  be  defined  in  all 

half  plan  t  >  0  because  that  |r1'(-«*,t)  _<  r  ,  |s^(-“,t)  |  <  s* 

+»  +°°  . 

by  Lemma  4.1  and  V  rft(*,t)  2Q.  ,  V  sn(*,t)  <  2Q.  by  (4.19). 

—00  —00 

By  Lemma  3.1  and  Lemma  3.2  we  can  obtain  the  existence  of  the  weak 
solution  (u,v) (x,t)  of  (0.1)  and  (0.2)  in  t  >  0  in  similar  way* 
to  that  in  theorem  3.1. 


mm*  - 
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Using  (4.9)  and  F(kh+0)  F(kh-O)  ,  k  -  1,2,...,  we  obtain 
F(nh+0)  <  (1- (1-c) oh) "f (0)  =  (1- (1-c) oh) nQL  , 


then  thus 
( 


+® 

u 

DV 

r  ( 

-  00 

+« 

DV 

sh( 

(4.29) 


It  is  not  difficult  by  the  means  similar  to  this  in  Lemma  4.2  to 
know  that 


+«° 


DV  rh(*  ,t)  <  DV  rh(* ,nh+0) 


,  for  nh  <  t  <  (n+l)h.  (4.30) 


+» 


DV  s1'(*,t)  <  DV  sh(*,nh+0) 


V 


let  i  tends  to  infinite,  by  (4.29)  we  have 
f 


DV 

—00 

r( 

+00 

DV 
—  00 

s( 

-(1-c)  at 

for  t  =  j2  and  k,j  are 


DV  s(*,t+0)<  Q1  e“(1_c)at 

•oo 

Therefore  by  (4.30)  and  (4.31)  we  obtain 

-  (1-c)  at 


positive  integers 


any 

(4.31) 


+«o 


DV  r(*  ,t)  <  Q,  e 
.00 


for  t  >  0 


(4.32) 


+» 


DV  <  Q.  e 


-(1-c)  at 


It  follows  from  (4.6)  that 


|rh(-®,t)-: 

|  sh(-~,t) 

therefore 

|r(-*»,t)-r 

|  s  (-00,  t)  -s 

Thus  then 
+® 

V  r(*,t)  < 
••00 

+00 

V  s(*,t)  < 
—00 

| r (x,t) -r |  < 
|s(x,t)-s  _< 

The  estimates 


i  <  | rQ (-«) -r  |  e”at  ,  |rh(+®,t) -r | <  | rQ (+®) -r| e”at 
I  1  |s0(-®)-s|e“at  ,  |sh(+«,t}-s|  _<  | Sq (+“) -s|e~at 

<  | rQ (-“) -r| e"at  ,  |r(+«,t)-r|  <  | rQ (+») -r | e”at 
|<  | sQ (-») -s| e~at  ,  |sh(+®,t)-s|  <  I sQ (+®) -s| e-at 


2Q1e-(1-c,ot 

+  (|r0(-®)-r|  + 

|rQ(+«)  t|  )e“at 

2Q1e'a-c)'xt 

+  (|sQ(-®)-s|  + 

ls0(+»)-s| )e“at 

+00 

v  r  ( • ,  t)  +  I 
—00 

rQ(— > -r | e"a 

+00 

V  s(*,t)  +  | 
—00 

I  sQ (-») -s [e~a 

• 

(4.25) , (4.26) , (4.27)  and  (4.28)  follow  from  (4.33) 
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5.  Entropy  Condition 


The  entropy  condition  for  quasilinear  hyperbolic  systems  of 
conservation  laws  has  been  extensively  studied  by  Lax,  P.D.  [3], 
Dafermos,  C.M  [1]  and  Liu,  Tai-Ping  [4].  Here  we  will  prove 
a  theorem  about  entropy  condition  for  the  weak  solution  of  (0.1) 
and  (0.2)  obtained  in  Section  3. 

Theorem  5.1.  Under  the  same  assumptions  of  theorem  3.1,  the 
weak  solution  (u,v) (x,t)  of  (0.1)  and  (0.2)  constructed  in 
theorem  3.1  satisfies  the  following  entropy  condition: 

||  fgt  E  +  gx  F  +  g(uU  -  i  V)>3xdt  >  0  (5.1) 

0<t<T 

for  any  non-negative  smooth  function  g  with  compact  support 

2 

in  0<  t  <T;  where  E(u,v)  =  ~  -  log  v,  F(u,v)  =  S. 

Proof.  Suppose  (u  .v  )  is  connected  to  (u0,v.)  by  a  1-shock  wave 

m  in  Xr  Jc 

with  the  speed  of  propagation  s^  on  the  left  in  the  solution  of 
(0.7),  let  G(V(vvt>)  =  ^(ECu^)  -  ElUj.Vj))  -  (FOi^)  - 
F(Ujl  vt^  '  ^hen  we  have 

G(V(W  -  0 
if  G('V(VV>  >  0 

XU 

because  s.  <  0,  v  >  v.  for  a  1-shock  wave  in  the  weak  solution 
x  in  x, 

of  (0.7). 


.u 


29. 


Therefore 


G(vm, (u£,v£) )  >  0  (5.2) 

Suppose  (u^,vr)  connected  to  by  a  2-shock  wave 

with  the  speed  of  propagation  s2  on  the  left  in  the  solution 
of  (0.7).  Let  G(vr,  )  »  s2(E  ur,vr)  -  E(um,vm))  - 

(F(ur,vr)  -  then  similarly  we  have 


"VW  >  0 


CS .  2  ■  1 


Using  (3. 3) (3. 4) (3.9)  and  (3.10),  we  have 


If  (gt  Eh  +  gx  Ph  +  g(uh  Uh  -  ^  ) dxdt  =  0(h)  +  3^  +  I2  (5.3) 


t<0<T 


Here 


Ij  ■  l  g(x,nh) (E(uh,rh) (x,nh-0)  -  E (un,vn(x ,nh+0) ) dx 
n=l 

l  f  9(sm  n  i  lEh]  “  IphJ)dt' 

=even  > «  m,n,l 


.h  h, 


m+n=even 
i-1,2 


(5.4, 


(5.5) 


where  the  speed  of  propagation  of  the  i-shock  wave  in 

(u,v)  (x, t (m+1)  ,n,  (uh,vh)  (m+an)  Jt,nh-0)  ,  (uh,vh)  ( (m+2+an)  A ,nh-0)  , 
t'-t-nh.  [•]  denotes  the  jump  of  the  quantity  in  the  bracket 
across  the  shock  wave.  Eh  =  E( (uh,vh) (x, t) ) ,  Fh  =  F ( (uh,vh) (x, t) ) . 

Using  (3.1)  (3.2),  we  have 


I2  -  0(h)  +  I'2 


where  I ' 2 


l 


m+n=even 
i-1,2 


/b 

J0*(s«,n,i(Bc5^,> 


(F(u,v))dt' 


(5.6) 


30. 


By  (5.1),  (5. 2),  we  have 


1'2  >  0 


(5.7) 


The  entropy  condition  (5.1)  follows  from  (5.3) (5.4)  (5.5) 


(5.6)  and  (5.7). 


Q.E.D. 


L  a-ll/Afci:1 
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